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ABSTRACT 

Context. Clustering in the luminosity of the afterglows of gamma-ray burst has been reported in the optical and X-ray. 

Aims. We investigate the possibility that a clustering in the luminosity of the afterglows of gamma-ray burst exists in the near infrared 

(J, H, K bands). We use observations of events from 1997 to the end of 2007. 

Methods. We correct the gamma-ray burst afterglow light curve for distance effects and time dilation, and rescale all light curves to a 
common distance of z = 1. We used only observations of signals emitted in the near infrared (in the burst frame). 
Results. We observe a clustering identical to the one observed in the optical and similar to the one observed in X-rays. We thus 
confirm previous work made in optical wavelengths and set a constraint on the total energy of the fireball. 
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1. Introduction 

Long Gamma-Ray Bursts (GRBs) are located at cosmologi- 
cal distances (Metzgeret al. 1997). These events are associated 
with supernovae (e.g. Hi orth etal. 1 120031: IStanek et al.1 l2003h 
and derive from the endpoint of the stellar evolution of mas- 
sive stars. It is thus tempting to hypothesize a standard candle 
behavior as for type la-supernovae. Indeed, several empirical 
coiTelations between observable quantities have been discovered 
(lAmati et al. l2002tlCTiirlanda et all2004 . trig gering the search 
for a standard candle based on GRBs (see e.g. iGhirlanda et alTI 
12006 ). This could have a very high impact on studies of the cos- 
molo gy in the redshift range 1-15 (for a review, see iMeszaros I 
I2OO6)." However, the use of GRBs for such studies requires the 
establishment of a robust indicator of their distance based on 
their intrinsic properties. This, in turn, implies that we must fully 
understand the properties of GRBs and their luminosity function. 

During the BeppoSAX era, we discovere d a clustering of 
the X -ray luminosities of GRB afterglows (Boer & Gendre ' 



200c ). This clustering was confirmed later (Gendre & Boer 



20051). and was also observed in the optica l (INardini et al. 120061: 



Kann et al. Il2006t iLiang & Zhang Il2006h . To date, GRBs clus- 



ter in three different groups according to their optical and X-ray 
luminosities; however, 10 % of nearby outliers are apparently 
underluminous and do not cluster (Gendre et al. 2008). Because 
of selection effects, this population of underluminous afterglows 
could be the tip of the iceberg of a more numerous undetected 
population. 

The physical reason for the clustering in luminosity is still 
unknown. It has been proposed that the microphysical parame- 
ters of the fireball (which fix th e afterglow propert ies) can take 
on only a few possible values ( Gendre et al. Il200"8l) . To confirm 
this hypothesis, we must investigate the luminosity clustering at 
longer wavelengths (compared to optical bands), so that the fire- 
ball parameters can be constrained. In this Letter, we present our 



findings obtained using infrared observations. They are interest- 
ing because : 

- infrared observations are numerous and can be used as a sta- 
tistically significant sample to estimate the positions of char- 
acteristic breaks within the afterglow broad-band spectrum, 
such as the injection frequency v„,; 

- the infrared sky is l ess sensitiv e to extinction compared to 
optical wavelengths jPei Ill992h . Because the extinction law 
(and thus the extinction correcti on) in distant galaxies is not 
well known jStratta et al. 11200 4^. this lowers the flux uncer- 
tainties due to the extinction correction. 

While our sample is not large enough to draw definitive con- 
clusions, in section|2]we show that there is evidence of luminos- 
ity clustering in the infrared. We will discuss this result in the 
fireball model framework in section |3] 



2. Data reduction and analysis 

2.1. Sample selection and definition 

We selected our sample from all long GRB afterglows observed 
from Earth in the near infrared (J, H, K bands) b etween 1997 
and th e first of January 2008. As pointed out by iKlotz et aT] 
(l2008h . some discrepan cies may occur be tween publications of 
data in GCN circulars ( lBarthelrriv1II998h and in refereed jour- 
nals. Because the latter are more rehable, we considered only 
data extracted from refereed journals and discarded data from 
other works. The advantage of near infrared (compa red to opti - 
cal wavelengths) is that extinction is less important dPei Ill992h . 
We thus reduced the sample to bursts in which the observed af- 
terglow emission was produced (in the burst frame) in the near 
infrared. This is easily done by computing the emission band 
using the observation band and the burst redshift. However this 
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Fig. 1. Light curves of our sample of afterglows, without distance correction. We plot them using the three bands we chose (see 
text for details) : the 8000 - 10000 A band (left), the 10000 - 13000 A band (middle), the 13000 - 16000 A band (right). Each burst 
is represented by a combination of different symbols and colors. See the electronic version for colors. 
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Fig. 2. Light curves corrected for distance and cosmological effects. We use the same symbols and colors as in Fig.[T] They are 
plotted at three different wavelengths : 9000 A (left), 1 1500 A (middle), and 14500 A (right). See the electronic version for colors. 



last constraint implies a selection bias based on the redshift (see 
discussion). This reduced sample is listed in Table[T] 

This sample of light curves was then correcte d for the extinc- 
tion du e to the Milky Way. We used the work of ISchlegel et al. I 
( 119981) for this purpose. Because the ext inction due to the host 
is more uncertain (see e.g. lLi et al. II2008I) we have not corrected 
the host extinction. 

Finally, the sample was distributed into three bands, based 
on the rest frame emission wavelength: emission within the 
8000 - 10000 A range, the 10000 - 13000 A range, and the 
13000 - 16000 A range. These bands do not correspond to any 
photometric filter used for observations. We defined them so that 
they are broad enough to include several afterglows within each 
band. We used three bands instead of a single large one because 
this reduces the flux interpolation needed to obtain a final value 
at a single wavelength (see next section). In the case of GRB 
030528 , we corrected th e flux for the presence of a bright host 
galaxy jRau et al. II2004I) . Note however that the host galaxy lu- 
minosity is not strongly constrained, thus the error bars on the 
corrected flux are large. Lastly, we removed from all data the 
initial part of the afterglow (see Sec. [3]). These light curves with- 
out further corrections are presented in Fig.[T] 

2.2. Distance corrections 

We used the method presented in lGendre & Boer I (|2005|) to take 
into account distance and cosmological effects. This method 
rescales all light curves to a common redshift of z = 1, and 



thus does not require knowledge of the value of Hq. We used 
a standard flat universe with Qa - 0.73 for the distance cor- 
rection. Since we already took into account the energy shift due 
to the cosmological effects by computing the emission band, we 
we do not need to rescale the observed frequency to a theoret- 
ical freguency_(as was done previously in the X-rays, and as 
e.g. iNardini et al. I l2006l did in the R band). We computed the 
monochromatic flux at the central frequency of each band, nor- 
malizing the fluxes by taking into account the spectral shape of 
the afterglow (using published spectral indexes or a value of 1 .0 
in the case of an unknown value). The corrected light curves are 
presented in Fig.|2] 

3. Results and discussion 

3.1. Resuits 

As one can clearly see in Fig. [21 we observe a hint of clustering 
of luminosity into two groups. From the sample of 19 events, 14 
belong to a group of bright events, 3 to a group of dim events, 
and 2 are isolated. In the following, we will use the notation of 
Gendre et al. (2008) : bright events in X-ray (x) or optical (o) 
are labeled with a /, and dim events are labeled with a //. A burst 
bright in the optical and dim in X-rays may thus be labeled as 
xll-ol. Outliers are labeled with a ///. We n ote that the infrare d 
groups are identical to the ones observed by Ka nn et al. I (l2006l) : 
Liang & Zhang (2006); Nardini et al. (2006) : we will thus not 
introduce the label ir but use the label o to discuss the infrared 
luminosity clustering. This reinforce the hypothesis that the clus- 
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Table 1. Sample of GRB afterglows used. We indicate for each 
burst its redshift, the bands observed at Earth and the references 



for the light curves used. 


GRB 


Redshift 


Observation band 


Reference 


GRB 970228 


0.70 


H,K 


a 


GRB 980703 


0.97 


H,K 


T 


GRB 990123 


1.60 


K' 


c 


GRB 991208 


0.71 


K 




GRB 991216 


1.02 


H,K 


e 


GRB 000418 


1.12 


K,Ks,K' 


7 


GRB 000911 


1.06 


H,Ks 


s 


GRB 010921 


0.45 


J,H,Ks 


T 


GRB 011121 


0.36 


J,H,K,Ks 


/ 


GRB 020405 


0.69 


H,K 


J 
~ 


GRB 021211 


1.01 


H,K 




GRB 030329 


0.17 


J,H 


7 


GRB 030528 


0.78 


H,K,Ks 


m 


GRB 031203 


0.11 


J,H,K 


n 


GRB 050126 


1.29 


Ks 





GRB 050408 


1.24 


K 


7 


GRB 060614 


0.125 


J 


T 


GRB 070125 


1.55 


Ks 


r 


GRB 070306 


1.50 


K 


s 


GRB 071010A 


0.98 


H,K 


t 



" lGalamaetal"l£l 997), Fruchter et al. ( 1999), Galama et al. (2000) 
'' IVreeswiik et al. I U99ft) . ,Castro-Tirado et al. . t2001.) . .Bloom et al. , 

( 119981) 

' LKulkami et al. ' (' 1999'),'Holland et al. ' ('2004|) 

[Sagar et al. (2000), Castro-Tirado et al. ( 200ll) 

^ [Halpern et al. (2000), Berger et al. (200n. lSagar et al. I ( l200dl) 
f [klose et al. (2000), Goros abel et al ■ (i2003l) 
s [Lazzati et al (2001), Mase tti et al. I d2005h 

* [Price et al. (2002a) 

' [Price et al. (2002b), Greiner et al.1 j2003h 
' [Price et al. ( 2003), S tratta et al. I ( l2005h 

* [Holland et al. (2004) 
'^loometal. (2004) 

lRauetal.i (2004), Rau et al. ' ('2 0M) 

" iGal-Yam et a l ( 200 4), Malesani et al. I ( |2004|) 
" iBerger et al I (120051) 

P iFoley et al. ' llQQ ^ 

" [Delia Valle et alTj ( l2006l) . ICobb et al. I ( 120061) 
'■ [Chandra et al. (260?) 
' [jaunsenetal. (.2008.) 
' ICovinoetal. I ( 120081) 



tering of infrared luminosity is not spurious. The flux difference 
between the two group is a; 3 magnitudes one day after the burst. 
We observe that GRB 031203, which does not cluster with the 
other bursts in the X-ray, seems also not to cluster in the near in- 
frared data (see Fig.|2] right pane l). GRB 060614, a nother nearby 
and peculiar burst (see e.g. Man gano et al. ||2007|) also does not 
cluster in luminosity. We discuss these results in the following 
sections. 

3.2. Selection effects 

As pointed out bv lGendre et al. I (l2008l) . the clu stering is only ap- 
parent in the afterglow data, and the value of lWillingale et al.1 
(l2007h can be used to discriminate between early and late light 
curves. However, most of the afterglows of our sample were ob- 
served before the launch of Swift, and the Tg value is unknown. 



To be conservative, and at the expense of the sample size, we re- 
moved the data taken before 30 000 s after the trigger Thi s value 
is larg er than most of the Ta values listed in Willingale et al. I 
(|2007|) . To avoid possible bias near this temporal cut-off, we will 
discuss in the following values obtained one day after the burst. 

The use of three bands instead of a single one is motivated by 
the reduction of uncertainties resulting from flux interpolations 
when correcting our sample : an error of +0. 1 in the spectral in- 
dex leads to less than 1.2% error in flux (compared to 5% in the 
case of a single large band). This has no incidence on the ob- 
served clustering, but reduces the sample size within each band. 

A last and most critical selection effect can be seen in Table 
[U : we selected only nearby events (compared to the mean Swift 
redshift). The bright events group ( ol) has a mean redshift of 1 . 1 . 
Most of our sample is composed of BeppoSAX events, it is thus 
not surprising that this low mean redshift is comparable to the 
BeppoSAX result. However, Swift has shown that the true red - 
shift distribution of GRBs was different dJakobsson et al. l2006l) . 
This redshift bias is due to our light curve selection. Removing 
all events emitted below the I band (but observed in infrared) 
implies that we remove all events above a redshift of 2.5. Thus, 
we cannot retain the majority of Swift bursts within our sam- 
ple. Howe ver, this should not a ffect our results strongly. As dis- 
cussed by iGendre et al. I (l2008l) . there are no strong differences 
in the afterglow properties of very distant and less distant bursts 
; this bias does not impact our results. However, a possible way 
to solve this bias is to add further observations within the far 
infrared bands. 

3.3. Expectations from tlie model 

The difiference between the fluxes of group / and // is ^ 3 magni- 
tudes one day after the burst; the optical and infrared groups are 
identical. This makes us confident that this clustering in lumi- 
nosity is not spurious despite our small sample. In the following, 
w e will use the results fo r the optical luminosity clustering given 
bv lNardini et al. I (l2006h . These authors reported a flux ratio be- 
tween group / and group // of ~ 26. Assuming no spectral vari- 
ation between optical and near infrared wavelengths, this trans- 
late into a difference in magnitude of 3.2, consistent with the 
observed separation. This implies that the clustering continues 
towards lower wavelengths than optical bands. 

Previous works have shown from spectral energy distribu- 
tions of GRB afterglows that optical an d near-infrared data can 
be fit using single power laws (see e.g. IStarling Il2008l for one 
example of such works), or in other terms, that the injection fre- 
quency v„, is lower compared to the infrared bands. We thus have 



Vm < ViR (1) 

with viR ~ 2x10''* Hz. In lGendre et al. I(l2008l) . we expressed 
several constraints on the fireball parameters. Using these con- 
straints and Equation[ri we can express a constraint on the total 
energy of the fireball (the constants C are listed in Table |2]l : 

E,, > CA:"-'vf^^'" (wind). ^ ' 

Clearly, the constraint indicated in Eq.|2]is not relevant in the 
case of normal GRBs : even in the case of a wind environment 
with a very low density (says A» - 10 '*) for a burst of the group 
xl-ol, we have 1.6 x 10"^ < £53. As one can clearly see in Eq. 
|2] in order to set a better constraint on this value, we need to 
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Table 2. Constants indicated in equations |2] and [3] for the three 
classes of afterglows in the wind and ISM. 



group ISM Wind 





Eq.|2] 


Eq.|3] 


Eq.|2] 


Eq.m 


xl-ol 


2.2 X 10** 


4.0 X 10" 


1.2 X 10' 


3.8x 10'- 


xll-ol 


1.3 X 10' 


2.2 X 10'" 


4.5 X 10' 


8.4 X 10'° 


xll-oll 


3.4 X 10« 


3.3 X 10' 


1.4 X 10' 


1.6 X 10'" 



observe at frequencies of the order of 10'^ Hz, but still above the 
injection frequency. 

3.4. Impact on the distance estimation 

In iGendre et al. I (l2008l) . we proposed a method to estimate the 
distance based on the afterglow flux observed one day after the 
burst in the 2-10 keV band. However, this method has a degener- 
acy, because of the existence of several groups. This degeneracy 
cannot be resolved by optical observations. To clear this prob- 
lem, we might have used infrared information. However, the fact 
that we observe exactly the same groups in optical and near in- 
frared makes it clear that we should use longer wavelengths; 

v^ = CE-^^"'\-'/'' Hz {ISM) ... 
v„ = CE;1°"a:'''Hz (wind). 

We express in Eq.[3]the position of v„, dSari et aT ll998h for 
the different groups (the constants C are listed in Table |2|l. As 
one can see, the position of this frequency is diff'erent for each 
group. This implies that at a frequency below the v,„ value of 
group xll-ol (the lowest one), we should observe three distinct 
groups. Observations at such frequencies could be performed by 
the IRAM observatory or the future ALMA observatory. 

4. Conclusions 

We investigated the clustering in luminosity of Gamma-Ray 
Burst afterglows using near infrared (J, H, K band) observations 
of events observed between 1997 and the end of 2007. We cor- 
rected the gamma-ray burst afterglow light curve for distance 
effects and time dilation, and rescaled all light curves to a com- 
mon distance of z = 1. We only used observations emitted in 
the near infrared in the burst frame. We observed a clustering 
of luminosities and a flux ratio identical to the ones observed in 
the optical. This implies that the injection frequency is located 
below the near infrared bands. We confirmed previous works on 
the clustering of afterglow light curves. 

Such a clustering has no clear theoretical explanation yet. In 
order to explain it, we need to set a strong constraint on the two 
remaining parameters of the fireball, namely the total energy and 
the density of the surrounding medium. For the latte r, radio ob- 
serva tions can be used to impose constraints (see e.g. lFrail et al. I 
12006 ). For the former, observations from about 100 GhZ to 10^^ 
Hz could help. However, this observational window has not been 
widely used for afterglow observations. It may be interesting to 
obtain several observations at microwave and millimeter wave- 
lengths that can be coupled with Swift and robotic telescopes to 
perform broadband spetro-temporal fitting of the afterglow. 

Finally, the constraints on the fireball parameters may be re- 
lated to the progenitor properties or to other conditions (such 
as the central engine properties) that would imply a clustering. 
In order to check the influence of the progenitor, one may add 
short GRBs (whose progenitor is thought to be a compact object 



binary) to the current study, that is based only on long GRBs 
(strongly suspected to be caused by massive stars). A clustering 
within the same groups would support a physical explanation 
not related to the progenitor but rather to the central engine of 
the burst. 
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